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We present our results concerning a systematical analysis of helioseismic implications on solar structure and 
neutrino production. We find Y ph = 0.238 - 0.259, R b /R® = 0.708 - 0.714 and p b = (0.185 - 0.199) gr/cm 3 . In 
the interval 0.2 < R/ Rq < 0.65, the quantity U = P/p is determined with and accuracy of ±5°/ 00 or better. 
At the solar center still one has remarkable accuracy, AU/U < 4%. We compare the predictions of recent solar 
models (standard and non-standard) with the helioseismic results. By constructing helioseismically constrained 
solar models, the central solar temperature is found to be T = 1.58 x 10 7 K with a conservatively estimated 
accuracy of 1.4%, so that the major unceratainty on neutrino fluxes is due to nuclear cross section and not to 
solar inputs. 



1. Introduction 

Helioseismology allows us to look into the deep 
interior of the Sun, probably more efficiently than 
with neutrinos. The highly precise measurements 
of frequencies and the tremendous number of 
measured lines enable us to extract the values of 
sound speed inside the sun with accuracy better 
than 1%. Recently it was demonstrated that a 
comparable accuracy can be obtained for the in- 
ner core of the Sun pj . 

The present work summarizes the results of our 
group in the last year concernig a systematic anal- 
ysis of helioseismic implications on solar structure 
and neutrino production. 

We quantitatively estimated the accuracy of 
solar structure properties as inferred from the 
measured frequencies through the so called inver- 
sion method. This analysis provided the base for 
quantitative tests of solar models. These tests are 
briefly presented here, see |l]-||] for more details. 

Concerning the organization of the paper, af- 



ter a short introduction to helioseismic data and 
their interpretation (sect. |J), we shall discuss the 
accuracy of solar properties as deduced from he- 
lioseismology (sect. ||). 

Next (sect. ||) we shall confront with helioseis- 
mology the predictions of Standard Solar Mod- 
els (SSMs). These tests are really a big success 
of recent SSMs, including elements diffusion, so 
that one can gain further confidence in the pre- 
dicted central solar temperature and T c and neu- 
trino fluxes <I>i. 

However, one can take a somehow different at- 
titude. The richness of helioseismic data can be 
used so as to supply information on some of the 
less certain inputs of solar models {e.g. solar opac- 
ity) . This is the approach we shall pursue to build 
Hcliosesimically Constrained Solar Models (HC- 
SMs), to be discussed in sect. ||. 

From the analysis of HCSMs it comes out that 
T c can be predicted with an accuracy of about one 
percent, a result obtained even if some of the most 
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controversial assumptions of SSM are released. 

This result contradicts recently advocated non- 
standard solar models, build so as to significantly 
decrease T c . Actually we show explicitely that 
these models are in conflict with helioseismic 
data, see sects. [| and [?]. 

Future prospects and applications of helioseis- 
mology are presented in sect. |[ 

2. Helioseismic data and their interpreta- 
tion 

Traditional and still most important helioseis- 
mic observables are frequencies of normal modes 
of solar oscillations. This is unlike in geophysi- 
cal seismic sounding where the primary data are 
travel times of seismic waves. A possibility of 
measuring the wave travel times in the Sun has 
been discovered not long ago [Q . The travel time 
data have been successfully applied to probing 
subsurface flows [§,1). The possibility of sound- 
ing short lived phenomena localized in the sub- 
surface layers is the advantage over the sounding 
based on the frequency data. Here, however, we 
are interested in sounding mean and solar struc- 
ture and therefore we will not discuss further this 
new method. 

The frequencies of solar oscillations are de- 
duced from time series of the intensity or radial 
velocity data. The latter are much more fre- 
quently used. Long-time continuous monitoring 
is essential for precise frequency determination. 
Four ground-based networks of automatic tele- 
scopes devoted to helioseismic observations are 
currently operated. The BISON and IRIS § 
networks measure radial velocity variations for 
the unresolved solar images. With this method, 
only modes of low spherical harmonic degrees 
(£< 3) are detectable. The GONG §] and TON 
networks use imaging instruments allowing to 
detect modes up to £ w 250 and £ « 700, respec- 
tively. There are three seismic instruments on the 
board of the SOHO spacecraft which started op- 
eration at the beginning of 1996. The GOLF Jn| 
instrument is similar to that used by IRIS. The 
LOI instrument fll2|] , which measures intensity 
variation, has a capability for detecting modes 
with £ < 8. The third is MDI, which an imag- 



ing instrument |l_3j. In the continuous mode of 
operation it detects modes up to £ = 250. 

Data for seismic probing from the imaging in- 
struments are usually provided in the form of cen- 
troid frequencies, vt n , and splitting coefficients, 
iktn, for the £,n multiplets, where n denotes ra- 
dial order of the mode. The coefficients describe 
the azimuthal dependence in oscillation frequen- 
cies according to the following relation 

Vlnm — VlnQ + ^ afcf„'Pfc(m), 
k=l 

where subscripts n and m denote radial and az- 
imuthal orders, respectively, and functions Vk(jn) 
are fc-order orthogonal polynomials. Data from 
the GONG network and from the MDI instru- 
ment provide values of v\ n $ and six coefficients 
for over 2000 (In) multiplets. For smaller num- 
ber of multiplets the coefficients up to k = 36 are 
available in the data sets from the MDI instru- 
ment. 

Values of a^'s reveal that the interior rotation 
rate is not very different from the surface rate and 
that the Sun is a very slowly rotating and almost 
perfectly spherical star. This fact significantly 
simplifies interpretation of the data. In particu- 
lar, the quantities i// n0 , a 2 k+i,in, and a 2 k,in are 
separate probes of, respectively, the radial struc- 
ture, the differential rotation and the aspheric- 
ity The possibility of sounding internal rotation 
is of great importance for understanding mecha- 
nism of the solar magnetic activity changing in 
the 11-year cycle. The data reveal that the dif- 
ferential rotation observed in the surface persists 
through the convective zone (some outer 30% in 
the radius). The transition to essentially uniform 
rotation takes place in the narrow zone, which 
has been the suggested site of the solar dynamo 
fl4|| . An important result in the context of this 
review is that there is no evidence for the rotation 
increase toward the center. In fact, according to 
some studies (e.g. p5[ ) there is a seismic evidence 
for lower rotation rate in the core than in outer 
layers. In any case the dynamical effect of cen- 
trifugal force is certainly negligible. There is also 
no evidence for a significant role of magnetic field 
in the Sun's interior. 

We now focus on the centroid frequencies and, 
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Figure 1. Frequencies of solar oscillations de- 
termined from 360-day measurements with the 
MDI/SOHO instrument. The error bars have 
been magnified by factor 1000. 



Figure 2. Relative differences between solar 
(shown in Fig. |l|) and frequencies calculated for 
standard solar model, "model S" of Ref. [30] . 



as an example, on the data from the MDI/SOHO 
instrument. The frequencies with the error bars 
are shown in Fig. |l|. The accuracy, which at 
the low frequencies is ~ 10~ 5 , visibly deteriorates 
at higher frequencies. The peaks in this part of 
the oscillation amplitude spectra are considerably 
broader, which is connected with the shorter life- 
time of the higher frequency modes. The total 
number of frequency data in the figure is 2047. 
The n — branch represents fundamental modes 
which are horizontally propagating waves obey- 
ing the same dispersion as surface waves in deep 
water. Their frequencies are fully determined by 
mean density in the Sun. The n > 1 branches 
represent p-modes which are standing acoustic 
(pressure) waves. These waves begin their down- 
ward propagation from the surface, where they 
are reflected, as nearly vertical rays. As a result 
of increasing adiabatic sound speed (a), the rays 
become gradually more oblique and, eventually, 
are reflected at the surface where the sound speed 
satisfies the condition a(R)/R = 2nv/y/£(£ + 1). 

Only a small fraction of p-modes from Fig. [j] 
penetrates the energy producing core. At v rj 3 
mHz and 1 = 1, the inner reflection surface oc- 
curs at R — 0.06i? Q , which is close to the maxi- 
mum the differential 8 B-neutrino flux of and well 
beneath the maximum of the differential pho- 
ton flux. However, at a similar frequency and 
I = 5 the reflection takes place at R = 0.18i?o, 
where the photon flux is already more than 92% 



of the total solar luminosity. Even for the low de- 
gree modes the frequencies are mainly sensitive 
to the structure of the outer layers. The modes 
which are very sensitive to the core structure are 
g-modes. Their frequencies, however, are below 
0.43 mHz and we may only hope that the solar 
g-modes will be some times detected. 

In Fig. H we show the relative differences in 
frequency between the Sun and one of the stan- 
dard solar models. Similar patterns are found 
for other models. The differences are small in the 
sense that there is no ambiguity in assigning the n 
value. However, they exceed the measurement er- 
rors frequently by two orders of magnitude. Thus, 
they are very significant and must be explained. 
In the case of the f-modes almost all the difference 
may be removed by an adjustment of the model 
radius |Q. Solar seismic radius is by 0.045% less 
than adopted in the model. For higher order p- 
modes the difference between model and solar fre- 
quencies rapidly increases with £. This immedi- 
ately suggests that most of the difference between 
the model and the Sun must be localized in out- 
ermost layers. Higher I implies a higher R at the 
lower reflecting surface and therefore a greater 
confinement to outer layers. 

We know that our treatment of the structure 
and oscillations in outer layers is inadequate. The 
problem is how to take into account effects of 
vigorous and nonadiabatic convection. We ex- 
pect the problem concerns the layers above R = 
O.99i?0, where p-modes propagate almost verti- 



cally and the local properties are ^-independent. 
This implies that the part of the frequency differ- 
ence may be modeled in the form F(y)/I, where 
/ is the calculated mode inertia assuming uniform 
normalization and F should be determined from 
the data. It is found that for modern standard so- 
lar models this part dominates in Sv but we are 
interested only in the small contribution arising 
in the rest of the interior. 

In the deeper layers the neglect of dynami- 
cal effects of convection is fully justified. Fur- 
thermore, the nonadiabatic effects in oscillations, 
which have not been taken into account in cal- 
culation of frequencies used in Fig. |2[ are cer- 
tainly negligible below R — O.99i?0. In view of 
the small values of 5v/i>, linearization about the 
reference standard model seems justified. Thus, 
we may use the variational principle for adiabatic 
stellar oscillations to connect the frequency dif- 
ference to the differences in structural functions 
between the Sun and its model. Our aim is to 
determine these differences. In general, we have 
to consider simultaneously two unknown func- 
tions of R. For one the choices are differences 
in density, p, pressure P, or any combination of 
them or their derivatives. All such functions are 
linked by linearized hydrostatic condition. We 
use here U = Pj p. For the other thermodynamic 
function the choices are the adiabatic exponent 
Y\ = (91nP/91np) ac i or squared adiabatic sound 
speed, a 2 = UT\. The problem is simplified if 
we make use the Ti(P, p, Y) relation, where Y is 
the fractional helium abundance, from the astro- 
physical equation of state data, (see e.g. Ref. 
|l6|]). Since in the chemically inhomogeneous in- 
terior the we may safely neglect 5Ti, the problem 
may be reduced to determination of the function 
SU(R) and the number SY p h~ the abundance of 
He in the layers above the base of the convective 
zone. In this way the basic equation for seismic 
probing of the internal structure becomes: 

,Sv, f .„ 5U .„ _ „. F{v) 



(1) 




where j = (n,£). The kernel Kj and numbers Jj, 
Ij are easily evaluated in terms of the eigenfunc- 
tions for adiabatic oscillations in the standard so- 
lar model. 
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Figure 3. Relative differences in U — p/p be- 
tween the Sun and the model inferred from the 
differences in v shown in Fig. |^. The symbols 
with error bars give the averaged values with the 
Gaussian-like localized kernels. The horizontal 
error bars give the full width at half-maxima of 
the kernels. Vertical error bars directly reflect the 
observational errors. The dashed line represents 
the result of inversion by means of the regularized 
least square method in which SU /U is searched in 
the form of a superposition of cubic splines. Near 
the center results obtained with this method are 
unreliable. The value of S Y ph = (4.0±0.1) x 1(T 3 
was determined simultaneously with SU/U. 



This equation may be applied to all p-modes. 
Thus, we have a large set of the integral equa- 
tions to determine the two functions and the 
Yph- Methods of determination of SU/U were 
described in details in Refs. jlTj and ]lg| ). There 
are two distinct approaches. The first is the reg- 
ularized least squares (RLS) method which con- 
sists in a discretization in terms of known function 
and a determination of parameters by means of 
minimization of the x~squared fit plus an integral 
term that smooths artificial oscillation. The other 
is the optimal averaging method which consists 
in seeking linear combination of individual ker- 
nels which are close to possibly narrow Gaussians 
centered at specified distances R. Corresponding 
combination of frequency differences yields the 
averaged values of SU/U weighted with these lo- 
calized kernels. Fig. || shows that the results ob- 
tained with these methods agree in a good agree- 
ment. 

Once we know SU(R) we may determine Sp(R), 
SP(R) and SMr making no additional assump- 
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tions. Hence we get seismic values of these three 
parametrs in the whole solar interior. At this 
stage we can determine the value R at the base 
of the convective, Rb, as well as the temperature 
T(R) for R > R b . What we cannot find is T(R) 
Y(R) profiles in the radiative interior. In order 
to disentangle the two last functions, which are of 
course most interesting in the context of the so- 
lar neutrino problem, we have to use data on the 
opacities and the nuclear reaction rates. This has 
been done (see Ref. |l7j for references to the orig- 
inal work) but we will not pursue this way here. 
Rather, we rely on directly inferred quantities as 
constraints on solar models. 

We see in Fig. || that the relative differ- 
ences in U are less than 5 x 10~ 3 in the whole 
[0.05 — 0.95]i? Q range of R. This agreement is 
very unlikely accidental and, thus, we regard it as 
a confirmation of the standard model of the solar 
evolution. In the outer part of the convective zone 
the agreement is worse but it may be easily im- 
proved by admissible modifications in the descrip- 
tion of the convective energy transfer. Beneath 
R = O.O5i?0 we just do not have enough informa- 
tion to make any statement. The small differences 
in the intermediate region may be eliminated by 
introducing ad hoc opacity modifications which 
are within the uncertainty of its calculation p9[ . 
However, the bump in 8U/U near the bottom of 
the convective envelope may have another cause. 
Richard et al. |^0| showed that this feature may 
be removed if one allows a weak rotation-induced 
mixing of elements below the bottom of the con- 
vective envelope. The effect was considered as a 
possible explanation of the deficiency of lithium 
in the Sun's atmosphere. Mixing brings lithium 
to deeper layers where it burns. It also results in 
a reduction efficiency of the gravitational settling 
leading to somewhat higher Y p h in the convective 
zone. Thus, the same effect may also explain the 
seismic correction to the photospheric He abun- 
dance {6 Y ph = 4 x 10~ 3 ). 

3. How accurate are solar properties as in- 
ferred from helioseismology? 

We performed a systematic and possibly ex- 
haustive investigation of the uncertainties of the 



helioseismic approach, in order to estimate the 
global error to be assigned to helioseismic deter- 
minations of solar properties. With this spirit, we 
analyse several physical quantities Q characteriz- 
ing the solar structure. Concerning the outer part 
of the sun, we discuss the photospheric helium 
abundance Y p h, the depth of the the convective 
envelope Rb, and the density at the bottom of the 
convective zone pb- Then we consider the "inter- 
mediate" solar interior (x—R/Rq = 0.2 — 0.65), 
analysing the behaviour of the squared isother- 
mal sound speed, U=P/p. Finally we investigate 
the inner region (x < 0.2), where nuclear energy 
and neutrinos are produced. 

We remind -see the previous section- that he- 
lioseismology measures only the frequencies {v} 
of solar p-modes, and quantities characterizing 
the solar structure are indirectly inferred from the 
{z/}'s, through an inversion method. Schemati- 
cally, the procedure is the following: 

a) One starts with a solar model, giving values 
Qmod and predicting a set {v mo( i\ of frequencies. 
These will be somehow different from the mea- 
sured frequencies, vq ± A^q 

b) One then searches for the corrections 8Q to 
the solar model which are needed in order to 
match the corresponding frequencies {v mo d + Sv} 
with the observed frequencies {^©}. Expression 
for b~v are derived by using perturbation theory, 
where the starting model is used as a zero-th or- 
der approximation. The correction factors SQ are 
then computed, assuming some regularity prop- 
erties, so that the problem is mathematically well 
defined and/or unphysical solutions are avoided. 

c) The "helioseismic value" Q@ is thus deter- 
mined by adding the starting value and the cor- 
rection []: 

Qq = Qmod + $Q . (2) 

For each quantity Q we have determined the 
partial errors corresponding to each uncertain- 
ties of the helioseismic method. In fact, there 
are three independent sources of errors in the in- 
version process: 

1 Concerning notation, we remark that 8Q indicates the 
correction to solar model to obtain helioscimic value (see 
Eq. whereas AQ indicates the estimated uncertainty 
on Q. 



G 



Table 1 

For the indicated quantities Q we present the helioseismic values Qq and the relative errors AQ/Q. All 
uncertainties are in °/ 00 ■ In the fifth and sixth row, for U = P/ p the values of the uncertainties are the 
maxima in the indicated interval. In the last two rows the results on U at points representative of the 
7 Be and 8 B neutrino production are shown. 
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i) Errors on the measured frequencies, which - 
for a given inversion procedure - propagate on 
the value of Qq . 

ii) Residual dependence on the starting model: 
the resulting Qq is slightly different if one starts 
with different solar models. This introduces an 
additional uncertainty, which can be evaluated by 
comparing the results of several calculations. 

iii) Uncertainty in the regularization proce- 
dure. Essentially this is a problem of extrapola- 
tion/parametrization. Different methods, equally 
acceptable in principle, yield (slightly) different 
values of Qq . 

It has to be remarked that, in view of the 
extreme precision of the measured frequencies, 
Avq/vq 10~ 4 pl|-p4]|, uncertainties corre- 
sponding to ii) and iii) are extremely important. 

For deriving a global uncertainty, we took a 
very conservative approach. May be that the pa- 
rameter variation was not exhaustive, and what 
we found as extrema are not really so, but actu- 
ally are quite acceptable values. In view of this, 
let us double the interval we found and interpret 
±(AQ)fe, as partial errors. Furthermore, let us 
be really conservative assuming that errors add 
up linearly. In conclusion, this gives: 

AQ = ±£|(AQ) fc |. (3) 

k 

In the following sections, we shall use this error 
estimate. 



3.1. Properties of the convective envelope 

Three independent physical properties of the 
convective envelope are determined most accu- 
rately by seismic observations, see Ref.(l): Y P h, 
Rb and p b |. 

The helioseismic predictions, together with 
their (conservatively estimated) accuracy, are 
shown in Table [l]. We remark that for all these 
quantities the uncertainty resulting from propa- 
gation of the frequency measurement errors is of 
minor importance with respect to the "system- 
atic" errors, intrinsic to the inversion method. 

3.2. The intermediate region 

The essential output of helioseismology is the 
reconstruction of the adiabatic sound speed pro- 
file, a. Our discussion is in terms of the related 
quantity U; as well known (see also sect. ||), 
a 2 =Y 1 U and Ti is given by the equation of state 
with an accuracy of 10 -3 or better (even the sim- 
plest EOS, fully ionized perfect classical gas yields 
Ti = 5/3 with an accuracy of about 10~ 3 ). 

By using the RLS method (see sect. ||) it is 
possible to derive directly the profile of U as a 
function of the radial coordinate throughout all 
the sun, except for the inner region (x < 0.1). 
The results (values of U and global errors), are 
summarized in Table [Ll 



2 A fourth seismic "observable", the sound speedat the 
convective radius is traditionally considered, e.g. J25| . We 
have not included it in our list since, as shown in Ref . [pj| , 
it is not an independent one. 
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Figure 4. The estimated global relative uncer- 
tainty on U — P/p (thick line) and that due to 
the observational errors (thin line). 



It is convenient to consider an intermediate so- 
lar region: 0.2 < x < 0.65. The upper limit is 
established by requiring that it is well below the 
transition to the convective zone, which we dis- 
cussed above. The lower limit is chosen so as to 
exclude the region of energy production (see next 
subsection). For this region the following com- 
ments are relevant: 

a) Each of the individual uncertainties nowhere 
exceeds 2 °/ 0o • 

b) Uncertainties from the accuracy on the mea- 
sured frequencies are of minor relevance with re- 
spect to the residual model dependence and to 
the sensitivity to the inversion parameters, see 
Fig. | 

c) All in all, even with the most conservative 
estimates, the helioseismic determination is ex- 
tremely accurate: \AU/U\ < 5°/ 00 throughout 
the explored region. 

3.3. The energy production region 

As well known, most of the energy and of so- 
lar neutrinos originate from the innermost part 
of the sun. According to SSM calculations, see 
e.g. Refs. about 94% of the solar lumi- 

nosity and 93% of the pp neutrinos are produced 
within x < 0.2, the region which we analyse in 
this section. 

Our results are summarized in Fig. |i| and Ta- 
ble! Clearly the precision worsens in this region, 
due to the fact that p-modes do not penetrate in 
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Figure 5. The difference between U as predicted 
by selected solar models, U mo d, and the helioseis- 
mic determination, U@, normalized to this latter. 
The dotted area corresponding to (4r). SUN24 
is the "model 0" of Ref. @; FRANEC96 is the 
"best" model with He and heavier elements dif- 
fusion of Ref. |p9|| ; BP95 is the model with metal 
and He diffusion of Ref. JCD is the "model 
S" of Ref. M. 



the solar core, and consequently the information 
one can extract from available experimental re- 
sults is limited, but still important. 

For example, at the production maxima of 7 Bc 
and 8 B neutrinos (xbc — 0.06 and xb — 0.04 
according to our best solar model with diffusion 
|2g| l ) the global accuracy is still a 2%, see Table 
[^Even at x = the accuracy is 3.5%. In conclu- 
sion, helio seismology provides significant insight 
even on the solar innermost core. 

4. Helioseismology and SSMs 

The comparison between the predictions of a 
few recent SSM calculations and helioseismic in- 
formation is shown in Figs. | and §. 

Concerning the (isothermal) sound speed pro- 
file, see Fig. ||, all models look generally good. 
Also SUN24, a model which neglects elemental 
diffusion, passes this test. 

The study of convective envelope is illuminat- 
ing, see Fig. |[ All models neglecting elemen- 
tal diffusion are in clear contradiction with helio- 
seismic constraint. On the other hand, calcula- 
tions where diffusion is included look in substan- 
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Figure 6. Hclioseismic determinations and so- 
lar model predictions of properties of the convec- 
tive envelope. The box defines the region allowed 
by helioseismology. Open circles denote models 
without diffusion, squares models with He diffu- 
sion, full circles models with He and heavier ele- 
ments diffusion, see Ref. 0. 



tial agreement with helioseismology. 

All this shows that the two approaches (profile 
of U and properties of the convective envelope) 
are complementary and both important. 

The previous arguments show that SSMs are in 
good shape. Actually, helioseismology provides a 
new perspective/definition of SSMs. 

Before the advent of helioseismology a SSM 
had three essentially free parameters, a, Yi„ and 
(Z/X)i n for producing three measured quantities: 
the present radius, luminosity and heavy element 
content of the photosphere. This may not look as 
a too big accomplishment, in itself. 

Nowadays, by using the same number of pa- 
rameters a SSM has to reprooduce many addi- 
tional data, such as Y p h, Rb, Pb, U(R), provided 
by helioseismology. 

Alternative solar models have to be confronted 
with these data too. 

5. Helioseismically constrained solar mod- 
els 

Helioseismic data are thus in agreement with 
recent SSM calculations, which use accurate 
equations of state, recent opacity tables and 
include helium and heavier elements diffusion 



HjlH-ijj, see also Ref. @. These SSMs 
yield central temperatures T c ^ssm which differ 
from each other by not more than 1%. However 
the uncertainties in the input parameters, mainly 
the opacity k and the heavy elements abundance 
Z/X, result in (AT C /T C ) SSM « 1 - 2%. 

From helioseismic observations one cannot de- 
termine directly temperatures of the solar inte- 
rior, as one cannot determine the temperature of 
a gas from the knowledge of the sound speed un- 
less the chemical composition is known. How- 
ever, it is possible to obtain the range of allowed 
values of the central temperature T c , by selecting 
those solar models which are consistent with seis- 
mic data. 

Our calculations are not model-independent, 
but we shall use in principle a wider class 
of models in comparison with SSMs, which 
we call helioseismically-constrained solar models 
(HCSM). These models are based on the same 
equilibrium and evolution equations as SSMs, but 
they differ in the choice of some input parame- 
ters. We generate this class of models by using 
the FRANEC code @ for the SSMs and varying 
the input parameters. Each choice of the set of 
parameters gives some value of T c . 

We obtain the range of allowed values of T c 
by selecting those solar models which are consis- 
tent with seismic data. More specifically, we shall 
determine the central temperature T c ^hcsm, as 
that of the model which gives the best fit to the 
seismic data and the uncertainties, AT c _hcsm, 
corresponding to the range spanned by models 
consistent with these data. 

We remind that the precise value of the tem- 
perature is governed essentially by two quantities: 
the radiative opacity k and the fraction of heavy 
elements Z /X. As well known the uncertainties on 
k and Z/X are of the order of 10%. Furthermore, 
these uncertainties do not correspond to clear ex- 
perimental or observational errors, rather they 
are determined by judicious comparison among 
published values. As an example, the uncertainty 
on k can only be estimated from the comparison 
among recent theoretical calculations. 

We allowed that both k and Z/X are rescaled 
by free multiplicative factors with respect to the 
value used in th SSM. These scaling factors are 
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Table 2 

Predictions of neutrino fluxes and signals in the CI and Ga detectors from HCSMs. Uncertainties 
corresponding to (AT /T)hcsm — ±1.4% are shown (first error) together with those from nuclear cross 
sections (second error). 

<F~^ [107cm7s] 4.81±0.53±0.59 

<f> B [10 6 /cm 2 /s] 5.96±1.49± 1.93 

CI [SNU] 8.4±1.9±2.2 

Ga [SNU] 133±11±8 



then determined by helioseismic constraints on 
the convective envelope. 

In this way (see Ref. |3^|), we obtained as the 
best estimate Thcsm = 1-58 x 10 7 K, with a con- 
servately estimated uncertainty (AT /T)hcsm = 
±1.4%, to be compared with the uncertainty of 
SSM calculations (AT/T) SSM = ±2.7%. 

We remark the following points: 

i) The "best" temperatures determined by means 
of helioseismology starting from different solar 
models converge. 

ii) The uncertainty (AT c )hcsm is (slightly) re- 
duced with respect to (AT c )ssm ■ 

iii) More important, the helioseismic uncertainty 
is related to observational data (whereas as noted 
above An, and thus (AT c )ssm, was derived just 
from comparison among theoretical calculations). 

Let us come over to neutrino fluxes, <&i (i—pp, 
7 Be, 8 B). Their dependence on the central tem- 
perature T c is parametrized as: 

$ 4 = ®i,SSM ( „ Tc ) • (4) 
\-<-c,SSAI / 

From numerical experiments with FRANEC, 
we found: (3 PP = -0.92, (3 Be = 7.9, B = 18. 

We have determined neutrino fluxes by starting 
with different SSMs, renormalizing their predic- 
tions to the same temperature Thcsm = 1.58 • 
10 7 K and to the same (updated) nuclear cross 
sections. The resulting fluxes and signals, all very 
close to each other, have been averaged to deter- 
mine the HCSM predictions shown in Table |[ 
where the first error corresponds to a conserva- 
tively estimate (AT c /T c )hcsm = 1-4%. and the 
scond one to a 3cr error on nuclear cross sections. 

We remark that nuclear physics uncertain- 
ties look larger than astrophysical ones. After 
the successful LUNA experiment measuring the 



3 He+ 3 He cross section |3^], the main uncertain- 
ties arise now from the astrophysical S-factor for 
3 He+ 4 He and p+ 7 Be reactions, which should be 
measured more accurately. 

6. A mixed solar core? 

Several authors have proposed non-standard 
solar models, where some ad-hoc mechanism is 
introduced, so as to reduce the central solar tem- 
perature by a few percent, which might alleviate 
(not solve) the solar neutrino puzzle. 

Clearly such a reduction of T c is in contrast 
with the conclusions of the previous section. Here 
we show explicitely that such models arc incon- 
sistent with helioseismic information. 

As an example, consider the case of a mixed so- 
lar core, firstly advanced by Shaviv and Salpeter 
in 1968 Hand recently revived by Haxton and 
Cumming |35|], see also Refs. ^8j,[56]]. Roughly 
speaking, mixing of H and He over an apprecia- 
ble portion of the sun enriches the innermost so- 
lar core with H. Hydrogen burning becomes more 
efficient and the solar luminosity is attained at 
lower core temperatures. This results in reduced 
T c , 7 Be+CNO and 8 B neutrino fluxes, again if 
a significant portion (R > O.li?©) of the Sun is 
mixed. 

Several different mixing processes (fast or slow, 
continuous or episodic) can be conceived |37| , j34| , 
H-0. While we refer to Ref. f§ for an extended 
discussion, we concentrate here on the case of fast 
continuos mixing. 

In Fig. we show (dashed line) the relative 
difference between the isothermal sound speed 
squared, U = P/p, as predicted by the mixed 
models and the helioseismic determination, ~Uq. 
The same quantity for our SSM is also shown 
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Figure 7. For the indicated values of R m ix, 
we present the relative difference between the 
isothermal sound speed as predicted by solar 
model with fast continuous mixing, U mo d, and the 
hclioscismic determination, Uq (dashed lines). 
The same quantity for our SSM is also shown 
(solid line). The dotted area corresponds to the 
uncertainty on Uq. 



(solid line). 

One sees a strong deviation of Mixed Core 
Models (MCMs), with respect to both Uq and 
Ussm, hi the mixing zone; this is a consequence 
of the change in "mean molecular weight" , /i, in- 
duced by mixing. In the approximation of perfect 
gas (accurate to the level of few per thousand in 
the solar core) one has U oc T/fi. Due to mixing, 
the innermost region is enriched with hydrogen, 
so that \i decreases (we observe that change of 
fx can be as high as 40%, whereas temperature 
change is at most a few per cent) and U increases. 
The opposite occurs near the edge of the mixed 
region. 

As the mixing area increases, the sound speed 
profile of the MCMs deviates more and more from 
the SSM prediction and it becomes in conflict 
with helioseismic constraint if R m i X > O.li?©. 

Neutrino fluxes predicted by MCMs are shown 
in Fig. ||. One sees that reduction of interme- 
diate ( 7 Be+CNO) and/or high ( 8 B) energy neu- 
trino fluxes is only achieved for Rmix > O.li?©. 

Similar conclusions hold for the other mixing 
mechanisms mentioned above, see Ref. ||. 

In brief, reduction of 7 Be and S B neutrino 
fluxes can only be obtained for extended mixing 
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Figure 8. The predictions of intermediate 
( 7 Be+CNO) and high energy ( 8 B) neutrino fluxes 
in solar models with continuous fast mixing, for 
the indicated values of Rmix- The prediction of 
our SSM (full diamond) is also shown. 



regions (R m ix > O.li?©), corresponding to solar 
models inconsistent with helioseismic constraints. 

We remind that the one percent accuracy of 
helioseismic information at R/Rq ps 0.1 is an es- 
sential ingredient for achieving this result. 

7. Helioseismology and physics of funda- 
mental interactions 

As another example illustrating the poten- 
tial of the helioseismic approach we discuss the 
p+p^d+e + +j/ e reaction. 

The rate of the initial reaction in the pp chain 
is too low to be directly measured in the labo- 
ratory (even in the solar center this rate is ex- 



tremely small, clearly of the order of 10 



yr 



- 1 ) 



and it can be determined only by using the the- 
ory of low energy weak interactions, together with 
the measured properties of the deuteron and of 
the proton-proton scattering. While we refer to 
Refs. Jl2],fl3],|7) for updated reviews, we remind 
that, as input of SSM calculations, one takes f43|j 
Spp,ssM = 3.89 • 10" 25 (1 ± 0.01) MeV b. 

We remark however that only theoretical esti- 
mates of S pp are available and observational in- 
formation would be welcome. In this respect, it 
is interesting to determine the range of S pp which 
is acceptable in comparison with helioseismology. 

In sect. |^ we pointed that there are two major 



uncertainties in building SSMs: solar opacity k 
and heavy element abundance £ =Z/X are only 
known with an accuracy of about ± 10%. By us- 
ing now k, C an d as S pp free parameters, we can 
determine the acceptable range of S pp by requir- 
ing that pb and Y p h are all predicted within 
the helioseismic range. The dependence of these 
quantities on k, ( and S pp has been determined 
numerically in Ref. p2| . 

Most of the information on S pp arises from data 
on p b as this observable depends strongly on S pp 
whereas it is weakly affected by the others pa- 
rameters. One can understand the dependence 
on S pp , at least qualitatively. A value of S pp 
larger than S PPi ssm implies smaller temperatures 
in the solar interior, which thus becomes more 
opaque (in other words, the region of partial ion- 
ization is deeper). Radiative transport therefore 
is less efficient and convection starts deeper in 
the Sun (Rt < Rb,ssAi), where density is higher 

(Pb > Pb,SSAl)- 

By restricting to the allowed ranges for the 
properties of the convective envelope, see Table 
, also taking into account the predictions of dif- 
ferent SSMs, we find: 

0.94 < S/Sssm < 1.18 (5) 

In conclusion, we remark that helioseismology 
provides the only observational constraint, al- 
though indirect, on the p+p — ► d + e + + v e re- 
action. 

Incidentally, we observe that this analysis com- 
pletely excludes the value S pp = 2.9S pp .ssm, re- 
cently found in Ref. H], by using a too rough 
description of the low energy pp scattering data, 
see also Ref. |pj. 

8. Future prospects and further applica- 
tions 

The helioseismic observation of the Sun will 
continue and the accuracy of the frequency data 
will improve. For solar physicists most hope 
is connected with the possibility of studying 
changes during the coming maximum of the mag- 
netic activity. Astronomers interesed in the an- 
gular evolution are looking forward to new data 
which should tell us whether indeed solar core 




R/R 

Figure 9. The estimated global relative uncer- 
tainty on p. 



rotates slower than the outer layers. This is so 
strange that most of them hesitate to believe it 
now. 

We do not expect revelations as far as inter- 
nal structure is concerned. The chances to find 
something that would shake the standard model 
of stellar evolution seem very small. Unquestion- 
ably, however, there is a room for improvement 
and helioseismology may help. 

As a few examples of relevance to calculation of 
precise solar models for evaluation of the neutrino 
fluxes as well to modeling evolution of other stars, 
we briefly outline a few topics we are presently 
working on. 

8.1. The solar density profile 

From the knowledge of U = P/p along the so- 
lar profile, by using the hydrostatic equilibrium 
equation, one can extract the solar density (see 
e.g. Jig]). We performed a systematic analysis of 
the uncertainties so as to provide a quantitative 
estimate of the accuracy A.p/p, the preliminary 
result being presented in Fig. ^. 

We remark that accurate predictions of the so- 
lar density, are important for determining the 
v c Vx transition probability in neutrino oscil- 
lation models where neutrino conversion is due to 
matter effects, see e.g. jlG] ff8| . 

8.2. Solar opacity 

As remarked above, solar opacity n(R) is the 
more uncertain ingredient of solar models. By 
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using the additional information provided by he- 
lioseismology one can try to transform it into an 
output of frequencies measurements. The idea 
is as follows: as we know U = P/p, we have 
one more information on solar structure, and we 
can use one of the equations of stellar equilib- 
rium/evolution for getting n. 

In fact, in the region beyond the energy pro- 
ducing core and below the convective envelope 
(0.2 < R/Rq < 0.65) one can use the radiative 
transport equation in the form: 



K(R)p(R) 



-lGnadT 4 
3i m ~d~R 



R 2 



(6) 



where L Q is the solar luminosity. In the same re- 
gion the perfect gas equation (for fully ionized H 
and He and for neutral Z) holds to high accuracy 
and the (almost) uniform abundance of He and 
Z can be inferred from the observed properties of 
the solar photosphera. In this way temperature 
gradients can be expressed in terms of gradients 
of U (again determined from helioseismology) so 
that the photon mean free path, A = l/(np), can 
be extracted, possibly with a few per cent accu- 
racy. 

This could provide an observational, although 
indirect, test of opacity calculation for plasma 
condition as in the solar interior, see also Ref. 
©■ 

8.3. The statistical distribution of nuclei in 
the solar core 

Already two decades ago, Kacharov et.al. |50t| 
and Clayton et.al. |5l[, speculated that the high 
energy tail of the distribution of protons in the 
solar core could depart from the Maxwell expo- 
nential form, 



dN Ma x{E) dN M 



dE 



dE 



-exp[-6(E/KT) 2 



(7) 



as an attempt to reduce the prediction of 8 B neu- 
trino flux by invoking a depletion in the number of 
protons with energy high enough to be captured 
by 7 Be nuclei. 

Recently interest in this matter was revived by 
Quarati et. al. |5i| , who provided some argument 
for such a depletion, in the frame of non extensive 
Tsallis statistics. 



In fact, sub barrier nuclear reaction provide a 
good laboratory for testing the high energy tail of 
the distribution, as they involve particles near the 
Gamow peak Eq, which is generally larger than 
the average thermal energy KT . Tiny deviation 
from the Maxwellian distribution, undetectable 
near KT, would be amplified at higher energies, 
see Eq. 0. 

If the energy distribution deviates from the 
maxwellian form, the reaction rates (av) would 
differ from those calculated with the Maxwell dis- 
tribution 

Clearly this also holds for the basic p+p— > 
e + +v e reaction. As we discussed in sect. |?j, its 
rate can be constrained by helioseismology. In the 
same way, possible deviations from the Maxwell 
distribution can be studied, see Ref. Ib3| 
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